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ABSTRACT 

We report the first results of a long term program aiming to provide accurate 
independent estimates of the Hubble constant [Hq) using the L(H/3) — a distance 
estimator for Giant extragalactic H ii regions (GEHR) and H ii galaxies. 

We have used VLT and Subaru high dispersion spectroscopic observations of a 
local sample of Hii galaxies, identified in the SDSS DR7 catalogue in order to re- 
define and improve the L(H/3) — a distance indicator and to determine the Hubble 
constant. To this end we utilized as local calibration or 'anchor' of this correlation, 
GEHR in nearby galaxies which have accurate distance measurements determined via 
primary indicators. Using our best sample of 69 nearby Hll galaxies and 23 GEHR 
in 9 galaxies we obtain Hq — 74. 3± 3.1 (statistical) ± 2.9 (systematic) km s~^Mpc~^, 
in excellent agreement with, and independently confirming, the most recent SNe la 
based results. 
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1 INTRODUCTION 

The accurate determination of the Hubble constant. Ho, is 
considered one of the most fundamental tasks in the inter- 
face between Astronomy and Cosmology. The importance of 
measuring the expansion rate of the Universe to high preci- 
sion stems from the fact that Ho, besides providing cosmic 
distances, is also a prerequisite for independent constraints 
on the mass-energy content of the Universe (e.g., Suyu et al. 
2012). 

The direct determination of the Hubble constant can 
only be obtained by measuring cosmic distances and map- 
ping the local expansion of the Universe, since the Hubble 
relation, cz — Hod, is valid and independent of the mass- 
energy content of the Universe only locally {z< 0.15). A va- 
riety of methods have been used to estimate -ffo, based on 
Cepheids, surface brightness fluctuations, masers, the tip of 
the red giant branch (TRGB), or type la supernovae [SNe la] 
(for general reviews see Jackson 2007; Tammann, Sandage 
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& Reindl 2008; Freedman & Madore 2010). In particular, 
the use of SNe la to measure the Hubble constant has a 
long history in astronomy (eg., Sandage & Tammann 1982; 
1990). The subsequent discovery of the correlation between 
the magnitude at peak brightness and the rate at which it 
declines thereafter (eg., Phillips 1993) allowed the reduction 
of the distance determination intrinsic scatter. However, one 
has to remember that SNe la are secondary indicators and 
their use relies on the determination of well-established local 
calibrators, like the Large Magellanic Cloud (LMC), Galac- 
tic Cepheids, the "maser" galaxy NGC 4258, etc. (cf. Riess 
et al. 2011). 

Indirect methods to measure -ffo have also been devel- 
oped (e.g. Bonamente et al. 2006; Suyu et al. 2010; Beutler 
et al. 2011), however, all of the indirect methods use as pri- 
ors other cosmological parameters, and thus the resulting 
Ho determinations are model dependent. 

Returning to the direct method to estimate Ho, an im- 
portant breakthrough occurred a decade or so ago by the 
HST Calibration program (Saha et al. 2001; Sandage et al. 
2006) who found Cepheids in local galaxies that host SNe la 
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and provided a Cepheid based zero-point calibration, and 
by the HST Key project (Freedman et al. 2001) who fur- 
nished a value of Hq = 72 ± 2(random) ± 7(systematic) km 
s~^ Mpc~^, based on Cepheid distances of external galax- 
ies and the LMC as the first rung of the distance ladder. 
This value was recently revised by the same authors, using 
a new Cepheid zero-point (Benedict et al. 2007) and the 
new SNe la of Hicken et al. (2009) , to a similar but less un- 
certain value of Ho = 73 ± 2(random) ± 4(systematic) km 
s~^ Mpc^^ (see Freedman & Madore 2010). Tammann et al. 
(2008) used a variety of local calibrators to recalibrate the 
SNe la and found a significantly lower value of Hq = 62.3 ±4 
km s^^Mpc^^. The difference has since been explained as 
being due to a variety of external causes among which the 
use of heavily reddened Galactic Cepheids and of less accu- 
rate photographic data (Riess et al. 2009a, b). 

The most recent analysis of Riess et al. (2011) uses new 
HST optical and infrared observations of 600 Cepheid vari- 
ables to determine the distance to eight galaxies hosting 
recent SNe la. The resulting best estimate for the Hubble 
constant is: Ho = 73.8±2.4 km s^^ Mpc^^ including random 
and systematic errors. 

From the above discussion it becomes clear that SNe la 
are the only tracers of the Hubble expansion utilized to-date, 
over a relatively wide redshift range (0<^ z<^ 1.5). Therefore, 
due to the great importance of direct determinations of the 
Hubble constant for cosmological studies (eg., Suyu et al 
2012) it is highly desirable to independently confirm the 
SNe la based Ho value by using an alternative tracer. 

H II galaxies have been proposed as such an alterna- 
tive. They are massive and compact (in many cases unre- 
solved) bursts of star formation in dwarf galaxies. The lu- 
minosity of H II galaxies is completely overpowered by that 
of the starburst. As a consequence they show the spectrum 
of a young Hll region, that indeed is what they are, hence 
their name. Their similarity with GEHR is underlined by 
the fact that the first examples of prototype H ll galaxies, 
I ZwlS and II Zw40, were called "Isolated Extragalactic H ll 
regions" and found to be observationally indistinguishable 
from GEHR in nearby galaxies (Sargent and Searle 1970). 
They are discovered mainly in spectroscopic surveys due to 
their strong narrow emission lines, i.e. very large equivalent 
widths. 

It is important to emphasise that the optical properties 
of H II galaxies are those of the young burst with almost no 
information (or contamination) from the parent galaxy. This 
is a direct consequence of selecting H ll galaxies as those sys- 
tems with the largest equivalent width (W) in their emission 
lines, i.e. W(H/3) > 50A. 

Because the starburst component can reach very high 
luminosity, Hll galaxies can be observed at large redshifts 
{z > 3). What makes these galaxies interesting cosmologi- 
cal distance probes (cf. Melnick, Terlevich & Terlevich 2000 ; 
Siegel et al. 2005) is the fact that as the mass of the starburst 
component increases, both the number of ionizing photons 
and the turbulent velocity of the gas, which is dominated 
by the star and gas gravitational potential, also increases. 
This induces a correlation between the luminosity of recom- 
bination lines, e.g. L(H/3) and the ionized gas velocity dis- 
persion a (see Terlevich & Melnick 1981; Hippelein 1986; 
Melnick, Terlevich & Moles 1988; Melnick, Terlevich & Ter- 



levich 2000; Fuentes-Masip et al. 2000; Telles et al. 2001, 
Bosch et al. 2002; Siegel et al. 2005; Bordalo & Telles 2011). 

A first attempt to estimate Ho, using Hll galaxies and 
GEHR as local calibrators, was presented in Melnick, Ter- 
levich & Moles (1988). The use of Hll galaxies as deep cos- 
mological tracers was discussed by Melnick, Terlevich & Ter- 
levich (2000) and Siegel et al. (2005) . Recently, we presented 
a thorough investigation of the viability of using H ll galaxies 
to constrain the dark energy equation of state, accounting 
also for the effects of gravitational lensing, which are ex- 
pected to be non-negligible for very high redshift 'standard 
candles' and we showed that indeed H ll galaxies can provide 
an important cosmological probe (Plionis et al. 2011). 

The aim of the current paper is to use Hll galaxies 
and a local calibration of the L(H/?) — o relation based on 
GEHR of nearby galaxies, as an alternative direct approach 
for estimating the Hubble constant over a redshift range of 
0.01 < 2 < 0.16. 



2 SAMPLE SELECTION AND OBSERVATIONS 

A sample of 128 H ll galaxies was selected from the SDSS 
DR7 spectroscopic data release (Abazajian, et al. 2009) 
within a redshift range 0.01 < z < 0.16, chosen for being 
compact (D < 5 arcsec), having large Balmer emission line 
fluxes and equivalent widths. A lower limit for the equiv- 
alent width of H/? of 50 A was chosen to avoid starbursts 
that are either evolved or contaminated by an underlying 
older stellar population component (cf. Melnick, Terlevich 
& Terlevich 2000). The redshift lower limit was chosen to 
minimize the effects of local peculiar motions relative to the 
Hubble flow and the upper limit to minimize any possible 
Malmquist bias and to avoid gross cosmological effects. 

In order to improve the parameters of the L(H/3) — a re- 
lation obtained from previous work, high-resolution echelle 
spectroscopy for the H ll galaxy sample was performed at 
8 meter class telescopes. We used the Ultraviolet and Vi- 
sual Echelle Spectrograph (UVES) (Dekker et al. 2000) at 
the European Southern Observatory (ESQ) Very Large Tele- 
scope (VLT) in Chile, and the High Dispersion Spectrograph 
(HDS) (Noguchi et al. 2002; Sato et al. 2002) at the National 
Astronomical Observatory of Japan (NAOJ) Subaru Tele- 
scope on Mauna Kea, Hawaii. The chosen setups provided 
UVES spectra centred at 5800 A with a slit- width of 2", 
giving a spectral resolution of ~ 22000. The HDS spectra 
were centred at ~ 5400 A, and with a slit width of 4" the 
spectral resolution obtained was ~ 9000. 

To obtain accurate total H/3 fluxes for the Hll galaxy 
sample, we performed long slit spectrophotometry at 2- 
meter class telescopes under photometric conditions and us- 
ing a slit width (8 arcsec) larger than the upper limit of 
the HII galaxies size in our sample. We used the Boiler & 
Chivens spectrographs at the 2.1 m telescope of the Obser- 
vatorio Astronomico Nacional (OAN) in San Pedro Martir 
and at the 2.1 m telescope of the Observatorio Astrofisico 
Guillermo Haro (OAGH) in Cananea, both in Mexico. 

Full details of the sample selection, observations and 
data reduction and analysis are given elsewhere (Chavez et 
al., in preparation). Here we summarize the relevant results 
regarding the determination of the distance estimator and 
Ho. 
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H^ and [O III] AA4959, 5007 line widths were measured 
fitting single gaussians to the line profiles. As previously 
found most H ll galaxies show line profiles that are well fit- 
ted by single gaussian (e.g. Melnick et al. 1988, Bordalo & 
Telles 2011). We cleaned the sample by first removing from 
the original list those Hll galaxies with either asymmetric 
or double/multiple line profile. We also removed those H ll 
galaxies showing rotation or large photometric errors in their 
H/3 fluxes or with an uncertain reddening correction. All this 
reduced the sample from 128 to 69 H ll galaxies. 

The values of the observed velocity dispersions, ao, were 
corrected for thermal (at) and instrumental (ai) broadening, 
and the final corrected dispersion was estimated according 
to: 



/ 2 2 



2^,1/2 



(1) 



The la uncertainties of the velocity dispersion were 
estimated from multiple observations computing the vari- 
ance over the repeated measurements; otherwise as the mean 
value of the obtained relative errors. 

H/3 integrated fiuxes were measured by fitting a single 
gaussian to the long slit spectra, while their la uncertainties 
were estimated from the expression (e.g. Tresse et al. 1999): 



ap = acD{2Np,^ + W/D) 



1/2 



(2) 



where ac is the mean standard deviation per pixel of the 
continuum on each side of the line, D is the spectral disper- 
sion, Npix is the number of pixels covered by the line and 
W is the line equivalent width. 

Heliocentric redshifts and their uncertainties were ob- 
tained from the SDSS DR7 and DR8 spectroscopic data 
when available, otherwise from our echelle data or the 
Spectrophotometric Catalog of Hll galaxies (Terlevich et 
al. 1991). The redshifts have been transformed from the 
heliocentric to the local group reference frames following 
Courteau & van den Bergh (1999) and corrected for the local 
bulk flow using the model of Basilakos and Plionis (1998). 
The la uncertainties were propagated using a Montecarlo 
procedure. 

To determine the zero point for the Z/(H/3) — a rela- 
tion, we obtained data from the literature for a sample of 23 
GEHR in 9 nearby galaxies whose distances have been mea- 
sured by means of well tested primary distance indicators. 

The details of the GEHR data will also be given in 
Chavez et al. (in preparation). For these objects, veloc- 
ity dispersions have been taken from Melnick et al. (1987), 
whereas distance moduli have been obtained averaging over 
the available measurements published after 1995, select- 
ing only those based on Cepheids, RR Lyrae, Mira vari- 
ables and eclipsing binaries except for those in IC 2574 and 
NGC 4236 for which only TRGB measurements are avail- 
able. The adopted distance moduli (/i) are listed as an inset 
in Figure[T] The global integrated H/3 fluxes and correspond- 
ing extinction were obtained from the values reported by 
Melnick et al. (1987). 



3 DETERMINATION OF Ho 

The procedure we use to estimate the Hubble constant com- 
prises three steps: 



log L(H p) 
r = 0.88, 
» NGC6822, M=23.36±0.01 
/ : M33, |U=24.49±0.01 
i NGC2366, M=27.68±0.20 
NGC2403, Ai=27.46±0.11 
IC2574, j4=27.90±0.02 
NGC4235, ;i=28.24±0.22 
MlOl, ;x=29.35±0.01 
- □LMC, ^=18.48±0.01 
SMC, ^=18.89±Q.01 



5.01 log a + (33.23± 0.05) 
^ 0.263, N = 23 




log cr(H/3) (km s"') 



Figure 1. L(H/3) — cr relation for the GEHR sample. The correla- 
tion parameters and the adopted individual distance moduli are 
given in the inset. The line is the best fit for the slope determined 
by the fit to the HII galaxies (see Eq. 3). 



log L(H/3) = ( 4.97±0.10)log a{Hf)+ (33.25±0.15) 
r = 0.96, rms = 0.236, N = 92 




1.2 1.4 

log cr(H^) (km s"') 



Figure 2. L(H/3) — ct relation for the joint Hll galaxies and GEHR 
samples. The least square fit considers the errors in both axes. 



(i) First we determine the slope of the I/(H/?) — a relation 
for Hll galaxies. Since the slope is independent of Hq we use 
an arbitary value of Hq to determine luminosities from the 
observed H/3 flux and the Hubble distancqj. 

(ii) We then determine the intercept of the relation from 
a flt to the 'anchor' GEHR sample, but flxing the slope to 
that determined in step one, i.e., that based on Hll galax- 
ies. Figure [T] shows the I/(H/3) — a relation for the GEHR 
sample. The slope of the correlation has been flxed to the 
value obtained from the Hll galaxies sample fltting in (i). 

The resulting L(H/3) — a correlation for the joint sample 
of GEHR and H ll galaxies is: 

logio L(H/3) = (4.97 ± 0.10) log^g a + (33.25 ± 0.15) (3) 



^ We have verified that the initial choice for the value of Hq does 
not alter the determined slope value. 
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Figure 3. Values of x for the grid of Hq. The solid line is a cubic 
fit to the points. The inset panel shows the value of x^ ~ Xmin ■ 



has r.m.s. logjQ L(H/3) — 0.236 and is shown in figure[2]. 

(iii) Finally we determine the value of Ho by minimizing, 
over a grid of Ho values, the function: 



x'{Ho[ 



= E 



[L,{a,)-U{Ho,f^,z,)f 
A?,,:+A?. 



(4) 



where the summation is over the Hll galaxies, at are the 
measured velocity dispersions (eq. [l]), Li{a) are the lumi- 
nosities estimated from the 'distance indicator' as defined 
in equation 3, Ala are their errors propagated from the un- 
certainties in a and the slope and intercept of the relation. 
Li{Ho, fi,Zi) are the luminosities obtained from the mea- 
sured fluxes and redshifts by using a particular value of Ho 
in the Hubble law to estimate distances, and Ai ^ are the er- 
rors in this last estimation of luminosities, propagated from 
the uncertainties in the fluxes and redshifts. 

Figure [3] shows the resulting x^ for the range of Ho values 
used, with the solid line being a cubic fit to the points. The 
la confidence limits of Hq were obtained from the values for 
which x^ ~ Xmin ~ 1 siucc the fit has only one degree of 
freedom (see the inset panel in Figure |3J. 

The value obtained for Ho using the above described 
procedure is: 



Ho = 74.3I3.J km s"^ Mpc" 



(5) 



Figure |?] shows the Hubble diagram for the sample of 
Hll galaxies used for the Ho value determination. The con- 
tinuous line shows the redshift run of the distance modulus, 
obtained from the linear Hubble law and the fitted Hq value, 
whereas the points correspond to the individual H ll galaxy 
distance moduli obtained through the L(H/?) — a correlation. 

The quoted Hubble constant uncertainty in equation 5 
reflects only the random errors, while systematic errors can 
also affect the mean value as well as the overall Hq uncer- 
tainty. We have identified as potential sources of system- 
atic errors the following: (a) the broadening of the emission 
lines, being contaminated by a rotational velocity compo- 
nent, (b) the internal structure/multiplicity of GEHR and 
Hll galaxies, (c) stellar winds affecting the line profiles (d) 
internal extinction, (e) coherent or peculiar motions affect- 



Figure 4. Hubble diagram for our sample of 69 H II galaxies. The 
thick points are the mean values for bins of 0.01 in redshift. The 
solid line shows the run with redshift of the distance modulus for 
Ho= 74.3. Residuals are plotted in the bottom panel and their 
distribution is shown in the inset. The r.m.s. value is 0.57 mag. 



ing the redshifts of nearby H ll galaxies, (f ) the age of GEHR 
and Hll galaxies, (g) the Malmquist bias, (h) variations in 
the IMF. The detailed discussion of these systematics is an 
important aspect of the Ho determination and will be pre- 
sented in a future paper (Chavez et al. in preparation). Here 
we briefly discuss these systematics, and the procedures used 
to minimize them. 

(a) and (b) To minimize the rotation and multiple com- 
ponent systematic effect we have used in our correlation only 
those objects with emission line profiles that are Gaussian 
and show no multiple components [for a discussion see Bosch 
et al. (2002)]. 

(c) The presence of weak extended (non-gaussian) wings 
in the emission line profiles introduces a small systematic 
effect. These weak wings are probably associated with stellar 
winds. The resulting effect is that taking into account the 
wings in the fit the final FWHM tends to be slightly smaller. 
This should affect similarly both GEHR and Hll galaxies. 
We estimate that this may introduce a systematic error of 
about 2 percent in Ho- 

(d) The extinction has been always estimated using the 
Balmer decrement method. We do not expect a sizeable sys- 
tematic effect associated with this correction. 

(e) H II galaxies tend to populate the voids so local pe- 
culiar motions should be relatively small. Furthermore, to 
minimize the effects of coherent bulk flows on the redshifts 
of the Hll galaxies, we imposed a lower radial velocity limit 
of 3000 km/s. In any case, we have also computed Ho in- 
cluding a local bulk flow correction and found no overall 
effect. 

(f ) The age of the GEHR and H ll galaxies affects their 
M/L ratio and therefore the zero point of the L(H/3) — a 
relationship. To minimize this effect we have speciflcally se- 
lected objects with W(H/3) > 50 A. This guarantees that 
the age of the star forming region is less than 6 Myr, thus 
minimizing the effect of evolution (Leitherer et al. 1999). 
We estimate that at most a plausible systematic difference 
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Table 1. Systematic error budget on the Hq determination 



Symbol 


Source 


Error (km s 


-1 Mpc-i) 


o-a,6 


Rotation, Multiplicity 


0.7 




CTc 


Stellar Winds 


1.1 




o-d 


Internal Extinction 


0.7 




'^f 


Object's Age 


1.4 




O-ff 


Malmquist Bias 


2.1 




o-h 


IMF 







Total 




2.9 





in ages between GEHR and H ll galaxies may affect Hq at a 
2% level. 

(g) We have calculated the Malmquist bias following 
the procedure proposed by Giraud (1987) adopting a power 
law luminosity function, with a slope a — —1.7. We have 
obtained a value of 2.1 km s~^ Mpc~^ aX z = 0.16, which 
we consider as one of the systematic error components. 

(h) The I/(H/3) — a distance estimator relies on the uni- 
versality of the IMF. Any systematic variation in the IMF 
will affect directly the M/L ratio and therefore the slope and 
zero point of the relation. The fact that our estimates of the 
Hubble constant are in agreement with those from SN la 
supports the hypothesis of a universal IMF. 

Table [T] shows the systematic error budget on the Ho 
determination. 



4 CONCLUSIONS 

It is indisputable that in the epoch of intense studies aimed 
at measuring the dark energy equation of state, it is of 
paramount importance to minimize the amount of priors 
needed to successfully complete such a task. One such prior 
is the Hubble constant Hq and its measurement at the ~ 1% 
accuracy level has been identified as a necessary prerequi- 
site for putting effective constraints on the dark energy, on 
neutrino physics and even on tests, at cosmological scales, 
of general relativity (see Suyu et al. 2012). Furthermore, it 
is highly desirable to have independent determinations of 
Hq, since this will help understand and control systematic 
effects that may affect individual methods and tracers of the 
Hubble expansion. 

It is within this latter strategy that our current work 
falls. We have carried out VLT and Subaru observations of 
a sample of nearby H ii galajdes identified in the SDSS DR7 
catalogue and 2m class telescopes spectrophotometry, in or- 
der to define their Z/(H/3) — g correlation, which we use to 
estimate the value of the Hubble constant. This is achieved 
by determining the zero-point of the distance indicator using 
GEHR in nearby galaxies, for which accurate independent 
distance measurements exist (based on Cepheids, RR Lyrae, 
TRGB and eclipsing binaries). 

Using our sample of 92 objects (69 Hll galaxies with 
z<, 0.16 and 23 GEHR in 9 galaxies with distances deter- 
mined via primary indicators) we obtain: 
J^o = 74.3±3.1(random) ±2.9 (systematic) km s"^ Mpc~\ 
in excellent agreement with, and independently confirming, 
the recent SNe la-based results of Riess et al. (2011). 
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